Two-band Eliashberg equations and the experimental T c of the diboride Mgi_ x Al x B 2 
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The variation of the superconducting critical temperature T c as a function of x in the diboride 
Mg\~ x Al x P>2 has been studied in the framework of the two-bands Eliashberg theory and tradi- 
tional phonon coupling mechanism. We have solved the two-bands Eliashberg equations using first- 
principle calculations or simple assumptions for the variation of the relevant physical quantities. We 
have found that the experimental T c curve can be explained only if the Coulomb pseudopotential 
changes with x by tuning the Fermi level toward the sigma band edge. In polycrystal samples the 
x dependence of the a and 7r-band gap has been found and is in agreement with experiments. 

PACS numbers: 74.20.Fg; 74.62.-c; 74.70.Dd 



The recent discovery of superconductivity at T=40 K 
in magnesium diboride Q] has stimulated intense inves- 
tigation, both from the theoretical and the experimental 
point of view. Now the electronic structure of MgB 2 is 
well understood and the Fermi surface consists of two 
three-dimensional sheets, from the tt bonding and anti- 
bonding bands, and two nearly cylindrical sheets from 
the two-dimensional a bands 2] . There is a large differ- 
ence in the electron-phonon coupling on different Fermi 
surface sheets and this fact leads to a multiband de- 
scription of superconductivity. The superconductivity in 
MgB 2 has been deeply studied in the past two years. 
More recently the substitutions of Mg with aluminum 
U have been investigated to understand the evolution 
of the pairing process by tuning the Fermi level toward 
the top of the sigma band edge. The Mg\_ x Al x B 2 al- 
loys show 0, H, l1 a continuous evolution through a 
complicated mixed phase from MgB 2 (x=0) to the fi- 
nal member AlMgB^ (x=0.5) where an ordered super- 
lattice structure of boron layers intercalated by alternat- 
ing layers of Al and Mg is formed. Even though the 
alloys with intermediate x are rather disordered, their T c 
is well defined and drops with decreasing x, becoming 
zero 0, 0, |(| for x > 0.5. We can see in Fig. 1 the ex- 
perimental data of critical temperature taken from ref. 
4-6. Two ranges of variation are present: a "high-TV 
range for < x < 0.3 that shows a slow variation with 
x, while a rapid variation occurs in the "low-T c " range 
0.3 < x < 0.5. This variation has been interpreted as to 
a 2D-3D cross-over, at x = 0.33, of the topology of the 
Fermi surface 0, . Theoretically, the transition tem- 
perature T c in Ai-substituted MgB 2 has been studied |7| 
within the two-band BCS formalism, as a function of the 
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FIG. 1: The value of the experimental T c (open circle) as a 
function of x, taken from ref. 4-6. 



aluminum content; the variation of the effective inter- 
band coupling as a function of x has been obtained from 
the experimental values of T c (x). Ab-initio calculations 
have been performed on ^-substituted MgB 2 by Profeta 
et aZ.0 and by de la Pena et aZ.0, within the virtual 
crystal approximation. These works correctly describe 
the experimental trends; in particular 8] , the behavior of 
the E2 g phonon frequency and of the critical temperature 
agree reasonably well with experiment. 

In this work we further investigate Mg\- x Al x B% (for 
< x < 0.5) by computing the variation of the crit- 
ical temperature as a function of x by means of the 
Migdal-Eliashberg two bands theory. For the peculiar 
characteristics of our system we must use a generaliza- 
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tion of the Eliashberg theory [lol ITU IT2| for systems with 
two bands 0, 0] . The two bands Eliashberg equations 
have alread y b een used for studying the MgB^ system 
13 0) ITsl jigj. Here we have ei ght input parameters: 
four (but only three independent |14|') electron phonon 
spectral functions oqJu))Fij(u>) and four (but only three 
independent) Coulomb pseudopotential ^jl*Au> c ) to repro- 
duce the experimental values of T c versus x by solving the 
Eliashberg equations. In this case there are four equa- 
tions to be solved for the calculation of the gaps A,(iw n ) 
and the renormalization functions Zi{iuj n ). In the ap- 
proximation of a flat normal densities of states and in- 
finite bands, the s-wave two bands Eliashberg equations 



are: 



uj n Zi (iuj n ) = oj n + irT 2J Aij (icj n - iu m )N j z (iu m ) + 
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FIG. 2: The spectral functions of the two-band model for the 
MgB2'. a a (solid line), nn (dashed line), an (dotted line) and 
7rcr (dashed dotted line), taken from ref. 15. 



Zi(iu! n )Ai(iuj n ) = 7rT^[Aij(io;„ - iu> m ) - /4y(w c )] ' 

m,j 

•6{\u! c \ — uj m )N^{iu) m 



where 6 is the Heaviside function, ui c is a cut-off energy, 
T lJ is the non magnetic impurity scattering rate in the 
Born approximation and 



Aij(iu) n - iu) m ) = 



50 dualjFicj) 

o (u„ - LU m ) 2 + UJ 2 



N J A (iu m ) = 



N z (iu m ) 



^LO^Z^iiuJra) + A 2 {iuj m ) Z 2 (iu) m ) 
^mZj(iu} m ) 

^jio 2 m Z 2 {iu m ) + A 2 (iuj m )Z 2 (iL0 m ) 



(3) 



(4) 



(5) 



where uj n — TtT(2n — 1) and n,m = 0, ±1, ±2.... In writ- 
ing the Eliashberg equations we have neglected the asym- 
metric part of the self-energy xd^n) (which is always 
equal to zero in the half-filling case) and the equation 
that represents the conservation of the particles neces- 
sary for calculating the shift of the chemical potential. 
These assumptions are correct if the normal density of 
states is symmetric or constant as function of energy and 
the bands are infinite. We have calculated the normal 
densities of state of the a and it band and, with good 
approximation, they are constant around the Fermi en- 
ergy for all x-values. Finally to use the infinite band 
approximation is correct for the 7r-band (half-bandwidth 
W„(x — 0) ~ 3 eV) while we can have some doubt for the 
cr-band {W a {x = 0) ~ 0.5 eV)H|. If we solve the Eliash- 
berg equations with these last finite values of the bands 



12, 24] for the MgB 2 (x = case), by using the same pa- 
rameters of the infinite bands approximation, we find T c 
=40 K. For having T c =39.4 K it must be ^(w c )=0.0333 
(see eq. 10) instead of 0.0326. Since in the first case T c is 
only slightly lower than the experimental one and /i(w c ), 
in the second case, is almost the same as in the infinite 
band case so we can conclude that to assume the bands 
as infinite is a plausible approximation. The situation 
can be different, of course, for large values of x. 

Finally, one can think that it would be simpler to use 
the Suhl formula |25| for explaining the T c versus x curve 
without solving the complete two-band Eliashberg equa- 
tions. The problem of the Suhl formula is the prefactor 
before the exponential: it is a phononic energy but in a 
two-band theory it is difficult to specify it exactly and so 
it is considered as another free parameter. 

For solving the two-band Eliashberg equations we must 
know how the six input quantities vary with x because, 
as shown in ref. 14, Ay (x)/Xji (x) — ^*j(x) / '^(x) = 

N 3 n (uj = E F ,x)/N l N (uj = E F ,x) where N J N (u = E F ,x) 
is the density of states at the Fermi level in the j-band 
and Ajj(x) are the electron-phonon coupling constants. 
For carrying out this task we must adopt some drastic 
approximations because almost all input quantities are, 
at the moment, difficult to determine. 

We begin with the four spectral functions of the MgB^ 
that we can see in Fig. 2 as calculated in ref. 15. We 
assume that the general shape of afjF(ui,x) does not 
change with x but only the Ay value, that is: 



Ay (x = 0) 



a 2 il F(uj,x = 0) 



(6) 



Since it is known that the details of a 2 F(u>) do not affect 
the resulting T c significantly |26|, this approximation is 
expected to be reasonable. 

From the definition of electron-phonon coupling con- 
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FIG. 3: Upper panel: calculated density of states at the Fermi 
energy Nn(u = Ef,x) in the <r-band (filled circles), in the 
7r-band (filled squares) and total (filled up triangles) as a func- 
tions of x. Lower panel: the phonon frequency uJE 2g obtained 
from experimental data as a function of x, see ref. 6. 



stant we have |2C 



N n (lu = E F ) < I 2 > 



Ami 



(7) 



M is the ion mass, the frequency flo is the frequency 
representative of phonon spectrum, Nn is the density of 
states at the Fermi level and < I 2 > is the average matrix 
element of the electron-ion interaction [2(J . The mass M 
is the boron mass [2lJ and does not depend on x. We as- 
sume that the average matrix element of the electron-ion 
interaction < I 2 > is, in the first approximation, constant 
because it is basically determined by the deformation po- 
tential which is almost || independent of x. We have cal- 
culated 8] , within the virtual crystal approximation, the 
dependence on x of N 3 N (see Fig. 3, upper panel) and 
we have identified the representative phonon frequency 
f^o with the E2 g phonon mode obtained by experimental 
data (see Fig. 3 lower panel). So we have that, the 
more important contribution to superconductivity in our 
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FIG. 4: Calculated electron phonon coupling constants Xij as 
a function of x in the four cases: aa (solid line), 7T7t (dashed 
line), an (dotted line) and na (dashed dotted line). 



system, the electron-phonon coupling constant X aa is: 

0) 



X aa (x) 



N^(u = E Fi x)ul(x 



N°(u; = E Fl x 



X acr (uj,x = 0) 



(8) 

In this way we assume that the change of Ei g in- 
fluence only the value of the electron-phonon cou- 
pling constant while, in the first approximation |2fi| |. 
we neglect its effects on the shape of the spectral 
function. For checking this hypothesis we have also 
used a^ .i 7 ' crcr (w, x) = L(u},Qq,T) — L(to,— £Iq,T) where 
L(lo, S7 , T) is a lorentzian curve with the peak in Qq(x) = 
loe 2 „( x ) an d half- width T(x) equal to the error bar of Fig. 
2 (lower panel). The result is almost the same as we can 
see in Figs. 5 (upper panel) and 6 (upper panel). 

For the other coupling constants, only for simplicity, 
we assume: 



\ij(x) 



Ef, x) 



N> N {ui=E F ,x = 0) 



Xij(x = 0) 



(9) 



with [HEil \ aa (x = 0) = 1.017, KA X = 0) = 0.448, 
X a7T (x = 0) = 0.213 and A 7r(T (a; = 0) = 0.155. 

In Fig. 4 we can see the calculated electron phonon 
coupling constants Ay as a function of x. We can see that, 
with our approximations, only \ aa depends strongly on 
x and, for x > 0.33, \ aa < A^. The Coulomb pseudopo- 
tential, calculated, for the first time, in ref. 15 and 16 is 

M 



^Lu c )N^ ot (E F ,p) 



2.23 



Merer M(T7r 
1 



N°(E F ,p) N°(E F ,p) 

1 2.48 
NZ(E F ,p) JV-(B F ,p) 



(10) 
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FIG. 5: Upper panel: the prefactor of the Coulomb pseu- 
dopotential /i as a function of x when the o? aa F{u, x) is a 
lorentzian (open circles), when is MgB2-like (filled circles) 
and when is MgB2-\ike and in presence of interband impu- 
rities, F """ — 1 meV, (open up triangles); in the insert the 
impurity scattering rate , necessary to reproduce the ex- 
perimental T c without change the MgB 2 physical parameters, 
versus x. Lower panel: the values of Coulomb pseudopoten- 
tial fi*j as a function of x in the four cases: a a (solid line), tvtt 
(dashed line), an (dotted line) and na (dashed dotted line) 
in the MgB 2 -\ike 

a%. a F(u},x) spectral function case. 

where fi(ui c ,x) is a free parameter and and Nffiiw = 
Ep, x) is the total normal density of states at the Fermi 
level. The numbers 2.23 and 2.48 in the Coulomb matrix 
have been calculated in the MgB 2 case and we suppose 
that, as first approximation, they do not depend on x 
i.e. the Coulomb pseudopotential depend by x only via 
the density of states at the Fermi level while an outer 
possible dependence by x, included the effect of the dis- 
order, is hidden in /i(w c , x) and is the same for the four 
values of fi*j (x) . In conclusion we put the cut-off energy 
u> c = 700 meV. So we have now only one free parameter, 
/i(w c , x) that must be determined by reproducing exactly 
the experimental critical temperatures. 

In our discussion we have neglected the effects of the 
disorder caused by Al doping because we don't know 
the dependence of r y by x and this fact lead to have 
a new free parameter in the problem. We seek of justify- 
ing this approximation and we examine two possibilities. 
One extreme case is one of thinking the Mg\- x Al x B 2 



FIG. 6: Upper panel: the value of Aj(iw n= o) of the <r-band 
(open circles) and of the 7r-band (open squares) calculated on 
imaginary axis, at T = T c /4, versus x when the (y^ a F(ui, x) 
is a lorentzian, when is M<;_B2-like (filled circles and squares) 
and when is MgB2-\ike and in presence of interband impu- 
rities, T aw = 1 meV (open up and down triangles). Lower 
panel: calculated rate 2Ai/fc.gT c for the cr-band (filled cir- 
cles) and of the 7r-band (filled squares) in the MgB2-like 
«ctct F(lo,x) spectral function case. 



as a contamination with Al impurities of the MgB 2 . 
In this case also we have only a free parameter for 
reproduce the experimental critical temperatures, T alT , 
because r CTCT and r™ don't affect |3 T c and T™ = 
N°(uj = E f )/N*(uj = E F )r a7r . In this way we find 
that T a7T growths fast in function of x as we can see in 
the insert of the upper panel of Fig. 5. The minimum 
critical temperature that we can obtain in this way is 
T c « 25.8 K and an isotropic gap A CT = A T « 4.1 meV 
but we cannot explain the lower critical temperatures. 
Moreover the authors of ref. 27 predict the T u > T lj w 1 
meV i.e. r y values found are completely to outside of 
the range expected. As second possibility we examine 
the case (2^ F a7T = 1 meV and we find the dependence 
of fx(u> c ) by x. As we can see in Figs. 5 and 6 (upper 
panel) the differences are almost imperceptible. For all 
these reasons we neglect the effects of impurities in the 
Eliashberg equations. 

In the case of MgB 2 (x = 0) we find fi(u> c ,x = 0) = 
0.0326 (uj c — 700 meV and maximum of energy equal to 
1 eV). In Fig. 5 we show the dependence on x of fi(u> c , x) 
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FIG. 7: The calculated rate of A^/A^ as a function of exper- 
imental T c . 



(upper panel) while in the lower panel we can see the 
values of n*j(x) calculated following eq. 10. We note a 
monotonic parabolic decrease until x — 0.4 while after 
it is present a rapid increase. This result shows that the 
coulomb pseudo-potential obtained in this way decreases, 
approaching the "shape resonance" determined by the 
superlattice of boron monolayers as predicted in several 
papers (HI- The reduction of the effective Coulomb re- 
pulsion in the pairing seems to be a driving term for 
rising the critical temperature. The final sharp increase 
can be due to an inadequateness of the model: probably, 
for x ps 0.5, the Migdal's theorem fails [2|| because there 
is a growth of the representative phonon frequencies and 
a strong reduction of the Fermi energy. 

After determining the dependence of the free parame- 
ter n{u c ) on x, we can calculate all the physical quantities 
by solving the Eliashberg equations. In Fig. 6 we show 
the value of A(iu> n —o), at T = T c /4, of the c-band and 
the 7r-band (solid circles and solid squares respectively) 
versus x. The values calculated are almost the same of 
those found by solution of the equation A? = A, (a; = A°) 
on the real axis with Pade approximants because we are 
in a weak coupling regime. In strong coupling regime 
instead the differences can be remarkable [3(J. We can 
see that the behaviour of the 7r-gap, as a function of x 
has a maximum for 0.164 < x < 0.198 while the c-gap 



decreases and, for x > 0.4, we find that A^ > A a . The 
experimental data HJexist only for T c > 20 K (x < 0.33) 
and are in good agreement with our theoretical predic- 
tions. For higher aluminum content, at present there are 
no measurements that confirm or deny the fact that the 
gaps merge for x ~ 0.4. In the lower panel we show the 
calculated ratio 2Ai/ksT c . Finally in Fig. 7 we can see 
the calculated rate of A n / A a as a function of experimen- 
tal T c (which is preferable due to the loch of agreement 
between the different experiments on the true aluminum 
content of the sample). 

It is worth to note that if, for the all four coupling con- 
stants (including the aa component), we use eq. 9, while 
the parameter /u(w c ) is fixed to the value (j,(u> c ,x = 0) = 
0.0326 used in the MgB 2 case, i.e. all input quantities of 
the Eliashberg equations depend by x only via the densi- 
ties of states at the Fermi level, we have a model without 
free parameters. This model produces a very good fit of 
the experimental critical temperatures for x < 0.33, but 
for larger x the predicted T c is much lower than the ex- 
perimental values. The value x = 0.33 corresponds a the 
transition from a pure two-dimensional Fermi surface to 
an overall three dimensional dispersion regime 8]. The 
behavior of the cr-gap is almost unchanged, while the one 
of the 7r-gap changes, and decreases monotonically. We 
have no arguments to justify the omission of the depen- 
dence of \ aa on uJE 2g - 

Further experimental work is necessary to determine 
the variation of the gaps with x in order to understand 
the correct theoretical approximation and to understand 
the variation of the pairing process tuning the Fermi level 
through the "shape resonance" condition. 

In conclusion we can affirm that not only the MgBi 
but also the same material doped with aluminum are 
weak coupling two band phononic systems well described 
by two-bands Eliashberg theory where the Coulomb pseu- 
dopotential and the intcrchanncl pairing mechanism are 
key terms to interpret the superconducting phase (it is 
the interband pairing that yields the superconductivity 
in the 7r-band). 

This research has been supported by INFM project 
PRA-UMBRA, INFM project RSPAIDSCS, INTAS 
project N. 01-0617 and CNR project "Progetto 5% Ap- 
plicazioni della superconduttivita ad alta T c " . 
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